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Abstract 


In this paper, the standard collocation approach is used to solve multi-order 
fractional integro-differential equations using Caputo sense. We obtain the 
integral form of the problem and transform it into a system of linear alge- 
braic equations using standard collocation points. The algebraic equations 
are then solved using the matrix inversion method. By substituting the 
algebraic equation solutions into the approximate solution, the numerical 
result is obtained. We establish the method’s uniqueness as well as the 
convergence of the method. Numerical examples show that the developed 
method is efficient in problem-solving and competes favorably with the 
existing method. 
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1 Introduction 


Fractional calculus is one of the subfields of mathematics that looks at the 
characteristics of the derivatives and integrals of noninteger orders. This 
discipline examines the notion and method of solving differential equations 
with fractional derivatives of unknown functions. In recent years, a signif- 
icant amount of interest in fractional calculus has emerged as a result of 
the fact that it may be used in a wide variety of fields of scientific inter- 
est; see [9]. Some of the numerical methods for the solution of fractional 
integro-differential equations developed in the literature include: Multi-order 
fractional by [12, 5, 16], Collocation method by [1, 3], Least square method 
by [13], Adomian decomposition method by [10], Chebyshev cardinal func- 
tions by [8], Laplace decomposition method by [11, 18, 14], Taylor expan- 
sion method by [7, 19], Haar wavelets by [4], Legendre Wavelets Method 
[6], and variational iteration method by [17]. Collocation approach to first- 
order Volterra integro-differential equations. The class of integro-differential 
equations was reformulated to assume an approximate solution in terms of 
the constructed polynomial. After solving for the unknown, we obtained a 
system of linear algebraic equations by collocating the resulting equation at 
various places within the range [0,1] [2]. The Laplace Adomian decompo- 
sition technique based on the Bernstein polynomial is employed to obtain 
an approximate solution for solving Volterra integral and integro-differential 
equations. Rani and Mistra [15] concluded that only orthogonal polynomi- 
als such as Legendre, Chebyshev, or Jacobi polynomials can improve the 
Adomian decomposition method. 


In this research, we present efficient method for solving multi-order frac- 
tional integro-differential equations with fractional derivatives of the form 


N b x 
DPy(a) = Dasa) ula) h(a) + | ker (0, t)y(t)dt-+ i ko(cr,t)y(t)dt 


(1) 


subject to the initial condition 
yD (aj) =r;, j=0,1,....n-1,nEN, Ban, (2) 


where y(a) is the unknown function to be determined, D® and D® are 
Caputo’s derivative, and h(x) is the force known prior. Moreover, k(,t) and 
k(a,t) are the Fredholm and Volterra integral kernel functions, respectively. 
Also, q;(x) is the known function and a; and \,; are known constants. 
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2 Basic definitions 


In this section, we present certain definitions and fundamental ideas of frac- 
tional calculus for the purpose of the formulation of the problem that has 
been presented. 


Definition 1. The Caputo derivative with order a > 0 of the given function 
f(x), x € (a,b) is defined as [11] 


1 
T'(m— a) 


Cey(z) = i “(e— 8)"-2-1y™ (s)ds, (3) 


wherem—1l<a<m, mE N, x >0. 


Definition 2. Let (a,),n > 0 be a sequence of real numbers. The power 
series in x with coefficients a, is an expression [11] 


N 
y(@) = ao + aya + age? +a30° +---ayaN = So ana” = g(a) A, (4) 
n=0 


where (x) = [1 a 22 ++» &%], A=[ap a, ++: ay)’. Then 
y(z,n) = a2"A, n=O0(1)N, ne Zr. 


Definition 3 (Standard Collocation Method (SCM)). This method is used 
to determine the desired collocation points within an interval, [a,b] and is 
given by [1] 


n=ap C—O ek eee A (5) 


Definition 4. Let y(x) be a continuous function. Then [3] 


Ny) 
olf (DE u(x) = (2) - Oar, (6) 


where m—-1<68<1. 
Definition 5. Let p(s) be an integrable function. Then [3] 


a “(n= s)P-2 S)das 
= ay [| @- 9" trea. (7) 


Definition 6. The Riemann—Liouville derivative of order a > 0 with n—-1 < 
a <n of the power function f(t) = t?~° is given by [11] 


T(p +1) 


pra 
T(p-—a-+l1) 


{Po (8) 


Definition 7. A metric on a set M is a function d: M x M —> R with the 
following properties, for all x,y € M [3], 
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(a) d(z,y) > 0, 

(b) d(z,y) =O => ar=y, 
(c) d(z,y) = dy, 2), 

(d) d(z,y) < d(x, z) +d(z,y). 


If d is a metric on M, then the pair (M, d) is called a metric space. 


Definition 8. Let (X,d) be a metric space. A mapping T : X —> X is Lip- 
schitzian if J a constant L > 0 such that d(Tx,Ty) < Ld(ax,y) for all x,y € 
X [3]. 


3 Mathematical background 


In this section, we develop an enhanced method for the numerical solution 
of multi-order fractional integro-differential equations. This method is based 
on the collocation approach and also considered power series polynomials as 
our basic function. 


Theorem 1 (Banach’s fixed point theorem). Let (X,d) be a complete metric 
space. It follows that each contraction mapping T’: X —> X has a unique 
fixed point x of T in X, such that T(x) = x 


Lemma 1 (Integral form). Let y(x) be a solution to (1) subject to (2). Then 
the integral form is 


« P(e) | [ogo tye Moat] as 
ie m5) [e- 5)e-2 ( [steno aa 
( 


4 Te fe — 3)8-1 [ kn(s,t)u(bd) ds, (9) 


where 


k=0 T(6 
Proof. Multiplying (1) by 918(-) gives 


of? (D°y(x)) = of 3 q(x 
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b 
+ olf (A(x) + off (/ (eta 


0 


+ gif ([ ka(s, bla) : (10) 


Using (6) on (9) gives 


kl L 
k=0 j=0 
b 
+ of2 (h(a) + off ( i, (esta 
0 
+ ol? ; 5 ‘ 
olf ( | ie u(t) (11) 


NY aOG 1 x = 
y(e) =>" a lee. ca (x — 5)? 


(8) Jo 
ee ee ae ; 
. a. ( ) ( ; ky( owe) d 
ay | ; , 
tray et Uf meow) 4 (12) 


k=0 

~ . M5 
x (Sota > 8 ya — (d(t)) us] ds 
+ _ [ (x — s)°"' h(s)ds + ne [ (2 — s)P? 


1 7 B-1 
2 ( in(esnotoa) dsA + aa, ie 8) 
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x ( | ; ka(s,t)0(t)t) dsA. (13) 


3.1 Method of solution 


Collocating at 2; in (13) gives 


N = 
ules) = Wes) + rom = aj) re) I er ae) 


[oe ore (600) it) ds A 
1 


dt™s 


Ly b 
pice acess re ts : 
ot. (ti — s) (/ Ki ( peta dsA 


a m5) i (a; — s)?-} ([ kn(s,t)o(t)tt) ds A, (14) 


where 


Ny 
= y (0) ky : x—s)?-'A(s S 
w( > il car. ( yo" h(s)ds. 
Simplifying (14) gives 
Thee Fim an ry 9 tal) 


” ie (s = ae ae! ae aan (d(t)) ir) ds 
= # 
TB) ies (x; = s)? 
x ( Je ki(s,t) ((t)) dt + f° ko(s,t) (6(0) dt) ds 
(15) 


b(xi)A = W(axi) + 


Factorizing the values of A from (15) gives 
yoy t fPee. —.9PNg, 
| $60) S30 TG a) rig 9 ai(s) 


” (6(t)) ar) dee 


x { fo(s— tyme} aaa 


1 price — g\6-2 
ria) fo 9) 


x ( Jo k(s,t) (o(t)) dt + f° ko(s,t) (@(t) dt) ds 


A=W(a;). (16) 
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Equation (16) can be in the form 


V(a,)A =W(2;), (17) 
where 
V (a) = (es) 5 eee [e — 8)?-*g;(s) 
D> OS 22 Tg — as) TB) Jo i 
i mj—aj—1 qm Ss — =a i t;—8s B-1 
(fem (6(0) at) ds — ae [ e-9) 
b s 
(/ ky (s, t) (oeyae+ [ ko(s,t) copa) ds (18) 
0 0 
and 
A=[a9 a: an] 


multiply both sides of (17) by V~!(a;) gives 


A =V~}(a;)W (2;). (19) 


Lemma 2. Let y(t) be approximated by (10) and let 


If gj(s) = s?/, then 


on) = PtP (n = a5 +0; +1) B+n—a; +P; 
tO agen ea ye +74 (21) 
Proof. Applying (3) and (7) into (20) gives 
N N Li 
18 | S> gj(x)D%y(2) | = st f(x — 8)? Ngy(s) 
; (3: ~~ > rm san A : 
[fe — yr-artyios) (at ds. (22) 
0 


Substituting (8) into (22) gives 


N 
of? (s a) 


j=0 
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N 
1 1 [ — 
= (x 3) SPi 
> Pim; — aj) T(8) Jo 


iz ie eae) ae ds A. (23) 


dt 
Let s—t = (1—v)s. Then t = vs 7 s dt = sdv. Substituting 
v 
them into (23) gives 
N 
of2 (37 aj(a)D%y(x) 
j=0 
N 
rT 1 1 sf 
> = | (a — s)®—1 gps 
ar IT(m; —a,;)P(n-—m,;+1)T(8) Jo 
1 
[srs | (1 —v)™— Vs dt | ds A. (24) 
0 


Simplifying (24), we get 


L(x:n) = P(n+i(n=a5+pjtl)gtn-astps , (25) 
T(n-—a,;+1P(6+n-—a; +p; +1) 


Lemma 3. Let y(t) be approximated by (10) and let 


b s 
E(x) = 918 / ka(s.du(dde + | baled : (26) 


If ky(s,t) = s"t? ko(s,t) = st”, then 


amie tC +4 1) Bir 
(ctn+ir@trtl 
Ren) = g A. 27 
(x; n) Tigtv+n+2) pbtgtvtn+] “a 
(utn+1)T(B+g+v+n+2) 


Proof. Applying (10) to (26) gives 
b s 
ff ki(s, t)y(t)dt ko(s, t)y(t)dt 
otf u(s.Qu(nar+ [ ko(s. ull | 


— a : rT—s B-1 ? : s : , 
7 ae. (m8) / ki ( du(oat+ | ko( swt ds: 
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Substituting ki(s,t) = s"t? ko(s,t) = s9t” gives 


b s 
of : ke (s,t)y(t)dt + | ite 


es “(a — )P} we 8 
=r | | ) (/ rata) 4 


+ ae [e — 8)? ( ste" y(ttt) ds. 


Applying (4) to (28) and simplifying give 


oI8 i" COO ee [ ko(s, t)y(t)dt 


0) 


1 x B-1 7 prre+l Ad 
=m | ee) |: ail 


ty fe 9 [e aa 
ae «Ss (rT S. 
T(8) Jo v+nt+1 


612 


(28) 


(29) 


Let « —s = (1—u)a. Then s = ux => ds = xdu. Substituting them into 


(29) gives 


ol? if ky(s, t)y(t)dt + ie ite 
0 
petetl 
_ T(8) (o+n+1 


T(8) (v+n+1) Jo 
Solving (30) gives 


betrtIT(r + 1) 


E(x;n) = air ae ala 


(vtnt+1)T(B+g+vu+n42) 


gh+r + 


Lemma 4. Let y(t) be approximated by (10) and let 


If h(s) = s™, then 
— Timi) Ba 
CO-Tetmty” 
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u) ajo" (ux)9t?tr+1 o du 


ghtgtutn+ 


A. (30) 


A. 


(31) 
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Proof. Applying (7) to (31) gives 


1 xz 
I (h -aa/ a — s)P-1s™ds 
ds 
Lett a—s=(l—u)a, s=uxr ras ds = xdu. Then 
U 


(32) 


Lemma 5. Let y(x) be the solution of (1) and (2). Then the numerical 
result gives 

y(x) = b(ai)V~*(ai) W (ai), (33) 
where 


Din+U0(n=ajt+pjtl)_p+n—aj+p; 
Tino, + UG +n =a) +p7+1)* 
prtet1p 1 
+ (r+ ) ght? 
(o+n+1)T(8+r+1) 


T(r FOTN 2) ghtrtotnt+l 
(otn+1)T(6+r+o4+n+2) 


V (ai) = 


and 


A agi) 
Wz.) = IE (0) «, V(m+ 1) pim 
a ak TB emt 


Proof. Approximate solution of (11) is 


From (19) we have A =V~!(x;) W(x;) where 


Din+UP(n=ajt+pjtl) pb +n—aj+p; 
T(n-a;+1)P(6+n-—a;+p;+1) * 
bret +1) 
'(otn+1)l(8+r41) 


i T(r TOrN+ 2) ghtrtetntl 
(otn+1)T(6+r+o4+n+2) 


ght 
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Substituting for A in the approximate solution gives the numerical result 


y(x) = b(@i)V~* (ai) W(2:). 


4 Uniqueness of the solution 


In this section, we establish the uniqueness of the method by introducing the 
following hypothesis: 


Ai: d= 
1:4 ee es 


Ay: ki = max |ky (a, t)| dt, 
«x€ [0,1] 


Hs :k3 = max i. |kin(a, t)| dt, 


x€[0,1] Jo 
A lw =") < Lm, lyn — yl, 
H. : = bm 
51U= mS = oaay 


Lemma 6. [q-contraction] Let T : X —> X be a mapping defined by Theo- 
rem | for y1, yg € X. Then T is q-contraction if and only if 


: ie! + Kk} +K3| <1 
T(8+1) [T(mj-a;t+1) t'? 
Moreover, there exist a unique solution of T. 


Proof. We have 


N x 
Tw)(2)=W0)+ > re aria f, ©) tal 


I= 


. “ll a OG bat ds + _ [oe =g)Pl 


b s 
x / n(s.Qn(bae+ [ itenntee ds 


= x ! ! ee s)?-19¢.(s 
(Tin) (= We) + oe any A | )ogs(s) 
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. [fc = yma (a ds + Re fe — s)P-1 


b s 
. / kn(s.t)yn(t)ae + | ital re 
0 0 


615 


(tm) @— Cw Ol=D ee aaty fe 9" Na 
«| femme fh — (0 at as 
0 


se 
tam fe-9 


fo lku(s.2)I ln) —ve()lat ] G, 
‘ Fe ee ae 


Taking maximum of both sides and using H; to Hs give 


(Tne) Tvee)) < ey |p + KE +] ay.) 


Since T' is a contraction, 


: | ua +Ki+K3] <1 
T(6+1) [T'(m; -— a; +1) 7 ? 


5 Convergence analysis 


In this section, we establish the convergence of the method by substituting 
the approximate solution into (3.0). We have 


1 1 
yn (2) = W(a) +>) T'(m,; — a3) (8) 


j=0 


« fe ay tay(0)| [ (o—m—ortafe? ae as 


= “(2 s)P72 , S S 
rand, | ) (/ ki ( punto) 


4: - [e _ 56-1 ([ ka(s,t)us (tt) de. (34) 
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Subtracting (9) from (34) gives 


Hence 


ce ae a eee 
<a f 9 Lae 


j — Oj) 


[eo - nme By (Oat ds 


[ow (s,t) Ew (t)dt| + 
0 


a ko (s,t) Ew(oal] ds. 


Therefore 


En (2) lho 
Ew (llc 


DP, ah 
<a [| ) 


The method of solution converges. 


ds. 


Dje0 Tim nayy 4 (8) [fo (s — t)"™- 09-4 de] 
a [I ky (s,t) dt + fy ke (s, t) dt] 


6 Numerical examples 


In this section, we present numerical examples to evaluate the effectiveness 
and clarity of the method. A MAPLE 18 program is used to perform the 
computations. Let y,(a) and y(x) be the approximate and exact solutions, 
respectively. Error y = |Yn(x) — y(«)|. 


Example 1. [6] Consider the following multi-order Fractional integro-differential 
equation: 


D¥Ty(e) = 22D! y(x) +2D°y(2)— ‘ee Dg@aee ‘ (a + 4) y(t)dt+ f(x) 


with this condition y' (0) = y(0) = 0 and exact solution y(a) = a? + #3, and 
_ (18) , re 5, (24 , 1 5 ©) .03 4.1. 

f(x) = (He tt) og? (1 Me) al ae 2° Taye 

xt a Tx 9 

12” 20. 12 20° 


Solution 1. Comparing with (1) and (2), 6 = 1.7,a, = 1.5,a2 = 0.5, ki (a, t) = 
(a +t), ko(x,t) = (a —-t). 


Using N = 3 for illustration. Applying (6) gives 
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1 1 * 1.7-1 ,2 
ule) = W(2) + TFs) aa a il 


[oe = pene et a dsA 


1 1 7 1.7-1 
* Fa—05) aa) oa a 


5 [fe S pros etal de A 


1 * 1.7-1 
fan. ee) 


( T'(n +1) ymin 4 I'(n+1) ae) 
T(n+m +1) T(n+m-+1) 
| Pet) mtn _ P(r +) «) dehy (90) 
T(n+m-+1) T(n+m+1) 
where 
= = : a — s)P—1 f(s)ds 
We) = Fa f @-9* "Yas. (36) 
Substituting f(s) into (36) gives 
1 . ae 
WO= Tay | (o~ aye? 
( T3) | FOR) )\ 2s, (FA@ _ TA) ) 3.5 
+ see a 5 
T5) | T@s5 T(2.5) | T@5) 
| TS) gO a ok | ae 
T(1.3) T(2.3) 12 20 12 20 


Simplify further 


_/7@) , T@)\ 1 oe 
ne (ri «ene gee 


(ro Spar ye 13-5 dg 
oe ra fe se oem 
ee 1 ;[ aye ds 


T(2.3) P(.7 
Lar can ae a 20P(1 a oe lle 


1 
( 
7 i 1.7-1 9 1 i, 1.7—1,,0 
my | cea ee 2011-7) J ae) a 
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(3) Pa5+1) 


_ (3) L742. 
va ee ": roa) Pirie = 
(BO, Te) POSH) nas 
T25) | 135) / 70743541) 
P@)_ TOS+1) arene 
TAs) P7034) 
TO. TOS) waa 
Tea Mir eis41) 
Di4+1) L744 P(5 +1) L7H 
nT 4441) Tr pad) a 
TL) ght 9 TO+)) L.7+0 (38) 
(iL? +141) 20F(1.7+0+1) 


Substituting (38) into (35) gives 


We obtain the result 


+0.9999999968x? + 1.0000000024x? 


_ fae x 10719 + 1.4273998650 x _ 


Table 1: Exact and approximate values of Example 1 


Exact N=3 


N=4 


N=6 


0.0781250000 0.0781249983 
0.3750000000 0.3749999996 
0.9843750000  0.9843749992 


0.0781249999 
0.3749999999 
0.9843749995 


0.0781250000 
0.3750000000 
0.9843749999 


2.000000000  1.9999999990 2.0000000000  2.0000000000 


Table 2: Absolute Error for Example 1 


x ERR3 ERR, ERRs (12]64 [6]7 
0.25 1.7e-9  1.0e-10 0.0 2,45e-4 —1.000e-5 
0.6 4.0e-10 1.0e-10 0.0 1.375e-3  1.200e-5 
0.75 8.0e-10 5.0e-10 1.0e-10 5.387e-3  2.13e-4 
1.0  1.0e-9 0.0 0.0 4.166e-3 8.970e-4 


Example 2. [5] Consider multi-order Fractional integro-differential equation 
of the form 
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619 
1 


1 
Dey(z) = -D'y(a)— y(e) + f y(at+e— 5 
0 
with the condition y(0) = y (0) = 1, and the exact solution is y(x) = # + 1. 
1 


Solution 2. Comparing with (1) and (2), 8 =2,~ 6=1.5,h(x) =a — 5. 


Use N = 8 for illustration. 
Write in the integral form 


ee) Fe 18) 5) [e 2 
[feo aanaee ra] ds A 
0 


T(n-2+1) 


coe “(a — 5)?1shdsA + ray |, @- 9" 


TQ 
| | ra ds A (39) 
0 
where , 7 ‘ 
W (a) = oo, (a—sy* (s _ 5) ds (40) 
ye) = We) + ae f w- 9 
T(n 2 1) n+0. 
Fe = 0.5)" ” | a 
ae 2-1.n Se i 2 —3)272 
- af @-9 s"dsA + roy I, ( ) 
ff ra ds A 
T(n+U0(n+1.5) 44 
We) We) +s ostmessy | 
T(n TT 1) n | T(n I 1) n 
Tin+3)- mad Tin+4)! ve oy) 
T(2) 5 AM). 4 (42) 


W@)= 5a" ~ 3B) 


for n = 0(1)N. Applying (41) and (42) gives 
y(2) = o(ai)V~* (ai) W (ai). 
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we obtain the result 


Ge 1.0000000000 + 1.0000000000a:-+ 
Y3\"7 * \ 88817841970 x 1071622 + 2.2204460493 x 10-1623 } ° 


Table 3: Exact and approximate values of Example 2 


x Exact N=3 Absolute Error 
0.2 1.2000000000 1.2000000000 0.00 
0.4 1.4000000000 1.4000000000 0.00 
0.6 1.6000000000 1.6000000000 0.00 
0.8 1.8000000000 1.8000000000 0.00 
1.0 2.0000000000 2.0000000000 0.00 


Example 3. [5] consider fractional fredholm integro-differential equations of 
the form 


1 
DMy(z) =D y(0) + [ ery(tjat + F(a), 
0 
where f(x) = e? — e**! with the condition y(0) = 0 and exact solution 
y(x) = e®. 


Solution 3. Comparing with (1) and (2), we have 8 = 1.5,a = 0.5,~ 
k(a,t) =e,~ f(z) =e" —e**!. 


Write in the integral form 


as 1 i es 
W0)=WO)+) meanest, & 3)! 


/ ee rma] ds A 
0 r 


m= mg +1) 
= neo [oe age? if e° rat ds A, (43) 
= Way= Teo ‘i "(2 — s)®-1 f(s)ds. (44) 


Use N = 8 for illustration. 
Substituting for 6 = 1.5,a = 0.5, f(a) = e® — e**! in (43) and (44) gives 


1 1 * 15-1 
ule) = We) + Bas) mas | ea 
. igeig Lae) 4 
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+ Tas [e- a [fe rat ds A 


W(a) = mas) i (g —s)'4-1 (e* _ err) ds 


for n = 0(1)N. Applying Lemma 4 gives 
T(n+1) 1 


y(x) = W(a) + Tin +05) (5) | (a — 5) 5-13 *15d5 A 


1 a peal Ss yn | 
+= u—s)” e° t"dt| ds A, 
T(1.5) a ) 0 


where 


Using (45) and (46) gives 
y(x) = $(ai)V~* (ai) W(ai). 


We obtain the result 


_ f 0.9990233401 + 1.0116982759x+ 
¥3 * \ 0.405067774922 + 0.3003042742x°3 


Table 4: Exact and approximate values of Example 3 


Collection-based numerical method ... 


(45) 


(46) 


Exact 


N=3 


N=5 


N=6 


1.2214027580 
1.4918246980 
1.8221188000 
2.2255409280 
2.7182818280 


1.2199681400 
1.4877329680 
1.8167324280 
2.2213811250 
2.7160936650 


1.2213960720 
1.4918178290 
1.8221150110 
2.2255348760 
2.7182828500 


1.2214031090 
1.4918249590 
1.8221190520 
2.2255409420 
2.7182817820 


Example 4. [16] consider the initial value problem of equation 


D?y(x) = «7D y(z) + 22D y(z) + 23 y(z)+ f(z), O<2<1, 
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Table 5: Absolute error for Example 3 


x ERR3 ERR; ERR¢« [5]N=18 

0.2 1.4346180e-3 6.686e-6 3.51e-7 0.21823e-7 
0.4 4.09173e-3 6.869e-6 2.6le-7 0.21586e-7 
0.6 5.386372e-3  3.789e-6 2.52e-7 0.86325e-7 
0.8 4159803e-3  6.052e-6 1.40e-8  0.12423-5 
1.0 2.1881634e-3 1.022e-6 4.60e-8 0.83792e-5 


f(x) = 6n'/? — 827/? = g/3 qi /2 


The exact solution is 
y(x) = ml 73 


Solution 4. Comparing with (1) and (2), we have 6 = 2,a; = 1.5,a2 = 
0.5, h(x) = 6rl/? — 87/2 — 18 y3 — g10/3q 1/2 


Use N = 8 for illustration. Applying (6) gives 


i 1 * 2-1 2 
ye) We) Seas) ny | ae al 


[eo = pas eet a ds A 


1 1 . 2—1 1/2 
D008) ay | (a — s)P* ab 


[eo - pros Teta ds A 


: . —s)?-! 51/3 [seq(s",n = 8 
+r [ & ) [seq(s”, 0,...,N)]ds A, (47) 


for n = 0(1)N, where 


1 e 16 
W(a2) = T(2) | (a—s)* (ons 8387/2 ; a iat) ds. 


Simplifying (47) gives 
a — Dn +P (n=1)0(3.5 +n) 
y(z) = W(@) + Fo 3)r(n 0.51 (5.5 +n) 


P(m=1) naa 4 , Pin +08) 
T(n +0.5) T(n + 2.8) 


6.5+n A 
gtt3.8 A, (48) 


6n'/?T (2) 81(3.5) ¢5 161(3) melas) ¢ 
WO-—To * tes)” ste) trey “© 
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Substituting (49) into (48) gives 


We obtain the result 


n= 


Collection-based numerical method ... 


—0.778452e — 42° — 0.212739e — we 


0.12350293e — 2x? + 1.76783815132° 


Table 6: Exact and approximate values of Example 4 


Exact 


N=3 


N=5 


N=7 


N=10 


0.0017725688 
0.0478593564 
0.2215710946 
0.6079910837 
1.2922026240 


0.0017002159 
0.0477585554 
0.2212000441 
0.6068809132 
1.2896573940 


0.0017710059 
0.0478450493 
0.2215300075 
0.6079068083 
1.2920556370 


0.0017729787 
0.0478592986 
0.2215712791 
0.6079913542 
1.2922025970 


0.0017721932 
0.0478593029 
0.2215710855 
0.6079910989 
1.2922026170 


Table 7: Absolute error for Example 4 


x ERR; ERR; ERR, ERR jo [16]10 
0.1 7.23529e-5 =1.5629e-6 4.099e-7 3.756e-7 7.45873e-7 
0.3 1.00801le-4 1.43072e-5 4.422e-6 5.35e-7 = -1.4833e-6 
0.5 3.7105e-4 4.1087e-5 1.845e-7 9.1le-9  1.74701e-6 
0.7 1.1101le-3 8.4275e-5 1.705e-7 1.52e-9 5.5116e-7 
0.9 2.54523e-3 1.46987e-4 2.7e-8 7.0e-9 2.47276e-6 


7 Discussion of results 


In this section, we discuss the numerical results obtained from the solved 
examples using the derived numerical method. 

In Example 1, the approximate solution obtained as N = 3 gives 
y3 = 1.8956614056 x 10~!° + 1.4273998650 x 10~12x + 0.9999999968a? + 
1.0000000024x?. Solving for N = 4 and N = 6, we obtained Table 1, which 
shows the results obtained from solving Example 1. Table 2 shows the ab- 
solute error of Example 1, and it indicates that as the values of N increase, 
the error becomes smaller and more consistent across all values of x. For 
instance, the least error of [12] at N = 64 is 2.45e — 4 while the least error 
in our method is 0.00 at N = 6. This confirmed that our method performed 
better. 
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In Example 2, the approximate solution obtained at N = 3 gives y3(a) = 
1.00000000 + 1.00000000002 + 8.8817841970 x 10~'6x? + 2.2204460493 x 
10-1°x3, which shows that the result converges to the exact solution as dis- 
played in Table 3. 

In Example 3, the approximate solution at N = 3 gives y3(x) = 
0.9990233401 + 1.0116982759x + 0.4050677749x? + 0.3003042742x3. Solv- 
ing N = 5 and 7, we obtained Table 4, which displays the results obtained at 
x = 0.2 to 1.0 for various values of N and the exact solution. The absolute 
error of Example 3 as shown in Table 5 indicates that as the values of N 
increase, the error becomes smaller. For instance, the least error in [5] at 
N = 18 is 0.21823e — 7 while the least error in our method at N = 6 is 
1.40e —8. This shows that the numerical method developed is consistent and 
converges faster. 

In Example 4, the approximate solution at N = 3 gives y3(x) = 
—2.5324187192 x 10~'8 + 6.5978333907 x 107122 — 1.0000000002x? 
+1.0000000000z?. Solving at N = 5, N =7, and N = 10, we obtained Table 
6, which shows the results obtained at x = 0.1 to 0.9 for various values of N 
and the exact solution. Table 7 shows the absolute error of problem 1, and it 
indicates that as the value of N increases, the error becomes smaller. We also 
compare our results with [16]. For instance, the least error in [16]Jat N = 10 
is 5.5116e — 7 while the least error in our method is 2.7e — 8 at N = 7. This 
clearly shows that our method performs better. 

Hence, from the numerical results obtained, we can conclude that the nu- 
merical method derived is efficient, consistent, and computationally reliable. 


8 Conclusion 


An enhanced numerical method was developed for the solution of multi- 
order fractional integro-differential equations with initial conditions using the 
collocation method. The numerical method derived is consistent, efficient, 
and reliable. Maple code was used to implement the developed method. 
Solved numerical examples showed that the method is reliable and suitable 
for such kinds of problems. 
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